EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 





CERN-PH-EP-201 2-003 
Submitted to: Eur. Phys. J. C 



(N 

o 

C 
(N 



Measurement of the Top Quark Mass with the Template Method 
in the tt — Hepton+jets Channel using ATLAS Data 



The ATLAS Collaboration 



X 

Oh 



> 

in 
in 
r^ 
in 

en 

o 

(N 

> 

X 



Abstract 

The top quark mass has been measured using the template method in the tt -> lepton+jets channel 
based on data recorded in 201 1 with the ATLAS detector at the LHC. The data were taken at a proton- 
proton centre-of-mass energy of y/s = 7 TeV and correspond to an integrated luminosity of 1.04 fb _1 . 
The analyses in the e+jets and ^+jets decay channels yield consistent results. The top quark mass is 
measured to be m top = 174.5 ± 0.6 stat ± 2.3 syst GeV. 
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Abstract The top quark mass has been measured us- 
ing the template method in the tt — > lepton+jets chan- 
nel based on data recorded in 2011 with the ATLAS 
detector at the LHC. The data were taken at a proton- 
proton centre-of-mass energy of ^/s = 7 TeV and cor- 
respond to an integrated luminosity of 1.04 fb _1 . The 
analyses in the e+jets and /i+jets decay channels yield 
consistent results. The top quark mass is measured to 
be m top = 174.5 ± 0.6 stat ± 2.3 syst GeV. 
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1 Introduction 



Template method 



The top quark mass (mtop) is a fundamental parameter 
of the Standard Model (SM) of particle physics. Due to 
its large mass, the top quark gives large contributions to 
electroweak radiative corrections. Together with preci- 
sion electroweak measurements, the top quark mass can 
be used to derive constraints on the masses of the as 
yet unobserved Higgs boson [1,2], and of heavy particles 
predicted by extensions of the SM. After the discovery 
of the top quark in 1995, much work has been devoted 
to the precise measurement of its mass. The present av- 
erage value of m t0 p = 173.2 ± 0.6 sta t ± 0.8 sys t GeV [3] 
is obtained from measurements at the Tevatron per- 
formed by CDF and D0 with Run I and Run II data cor- 
responding to integrated luminosities of up to 5.8 fb _1 . 
At the LHC, m t0 p has been measured by CMS in tt 
events in which both W bosons from the top quark 
decays themselves decay into a charged lepton and a 
neutrino [4]. 
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The main methodology used to determine m t0 p at 
hadron colliders consists of measuring the invariant 
mass of the decay products of the top quark candi- 
dates and deducing m t0 p using sophisticated analy- 
sis methods. The most precise measurements of this 
type use the tt — > lepton+jets channel, i.e. the decay 
tt —¥ £i>be <7iq2&had with £ = e,/x, where one of the W 
bosons from the tt decay decays into a charged lepton 
and a neutrino and the other into a pair of quarks, and 
where bi (&had) denotes the 6-quark associated to the 
leptonic (hadronic) W boson decay. In this paper these 
tt decay channels are referred to as e+jets and /z+jets 
channels. 

In the template method, simulated distributions 
are constructed for a chosen quantity sensitive to the 
physics observable under study, using a number of dis- 
crete values of that observable. These templates are fit- 
ted to functions that interpolate between different input 
values of the physics observable, fixing all other param- 
eters of the functions. In the final step a likelihood fit 
to the observed data distribution is used to obtain the 
value for the physics observable that best describes the 
data. In this procedure, the experimental distributions 
are constructed such that they are unbiased estimators 
of the physics observable used as an input parameter 
in the signal Monte Carlo samples. Consequently, the 
top quark mass determined this way from data corre- 
sponds to the mass definition used in the Monte Carlo. 
It is expected [5] that the difference between this mass 
definition and the pole mass is of order 1 GeV. 

The precision of the measurement of mt op is lim- 
ited mainly by the systematic uncertainty from a few 
sources. In this paper two different estimators for m top 
are developed, which have only a small statistical corre- 
lation and use different strategies to reduce the impact 
of these sources on the final uncertainty. This choice 
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translates into different sensitivities to the uncertainty 
sources for the two estimators. The first implementation 
of the template method is a one-dimensional template 
analysis (Id-analysis), which is based on the observable 
i?32 , defined as the per event ratio of the reconstructed 
invariant masses of the top quark and the W boson re- 
constructed from three and two jets respectively. For 
each event, an event likelihood is used to select the 
jet triplet assigned to the hadronic decays of the top 
quark and the W boson amongst the jets present in the 
event. The second implementation is a two-dimensional 
template analysis (2d-analysis), which simultaneously 
determines m top and a global jet energy scale factor 
(JSF) from the reconstructed invariant masses of the 
top quark and the W boson. This method utilises a x 2 
fit that constrains the reconstructed invariant mass of 
the W boson candidate to the world-average W boson 
mass measurement [6]. 

The paper is organised as follows: details of the 
ATLAS detector are given in Section 2, the data and 
Monte Carlo simulation samples are described in Sec- 
tion 3. The common part of the event selections is given 
in Section 4, followed by analysis-specific requirements 
detailed in Section 5. The specific details of the two 
analyses are explained in Section 6 and Section 7. The 
measurement of mtop is given in Section 8, where the 
evaluation of the systematic uncertainties is discussed 
in Section 8.1, and the individual results and their com- 
bination are reported in Section 8.2. Finally, the sum- 
mary and conclusions are given in Section 9. 



2 The ATLAS detector 

The ATLAS detector [7] at the LHC covers nearly the 
entire solid angle around the collision point 1 . It consists 
of an inner tracking detector surrounded by a thin su- 
perconducting solenoid, electromagnetic and hadronic 
calorimeters, and an external muon spectrometer in- 
corporating three large superconducting toroid magnet 
assemblies. 

The inner-detector system is immersed in a 2T axial 
magnetic field and provides charged particle tracking in 
the range \rj\ < 2.5. The high-granularity silicon pixel 
detector covers the vertex region and provides typically 

1 ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the 
y axis points upward. Cylindrical coordinates (r, <j>) are used 
in the transverse plane, <j> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 8 as r/ = — lntan(0/2). Transverse momentum 
and energy are defined as pt = psin# and Et = EsinO, 
respectively. 



three measurements per track, followed by the silicon 
microstrip tracker which provides four measurements 
from eight strip layers. These silicon detectors are com- 
plemented by the transition radiation tracker, which en- 
ables extended track reconstruction up to \rj\ — 2.0. In 
giving typically more than 30 straw-tube measurements 
per track, the transition radiation tracker improves the 
inner detector momentum resolution, and also provides 
electron identification information. 

The calorimeter system covers the pseudorapidity 
range \r]\ < 4.9. Within the region |?/| < 3.2, elec- 
tromagnetic calorimetry is provided by barrel and end 
cap lead/liquid argon (LAr) electromagnetic calorime- 
ters, with an additional thin LAr presampler covering 
\f]\ < 1.8 to correct for energy loss in material upstream 
of the calorimeters. Hadronic calorimetry is provided by 
the steel/scintillating-tile calorimeter, segmented into 
three barrel structures within \rj\ < 1.7, and two cop- 
per/LAr hadronic endcap calorimeters. The solid an- 
gle coverage is completed with forward copper/LAr and 
tungsten/LAr calorimeter modules optimised for elec- 
tromagnetic and hadronic measurements respectively. 

The muon spectrometer comprises separate trigger 
and high-precision tracking chambers measuring the de- 
flection of muons in a magnetic field with a bending 
integral up to 8 Tm in the central region, generated by 
three superconducting air-core toroids. The precision 
chamber system covers the region \t)\ < 2.7 with three 
layers of monitored drift tubes, complemented by cath- 
ode strip chambers in the forward region. The muon 
trigger system covers the range \r]\ < 2.4 with resistive 
plate chambers in the barrel, and thin gap chambers in 
the endcap regions. 

A three-level trigger system is used. The first level 
trigger is implemented in hardware and uses a subset 
of detector information to reduce the event rate to a 
design value of at most 75 kHz. This is followed by two 
software-based trigger levels, which together reduce the 
event rate to about 300 Hz. 

3 Data and Monte Carlo samples 

In this paper, data from LHC proton-proton collisions 
are used, collected at a centre-of-mass energy of ^/s = 
7 TeV with the ATLAS detector during March-June 
2011. An integrated luminosity of 1.04 fb _1 is included. 
Simulated tt events and single top quark produc- 
tion are both generated using the Next-to-Leading Or- 
der (NLO) Monte Carlo program MCONLO [8,9] with 
the NLO parton density function set CTEQ6.6 [10]. 
Parton showering and underlying event (i.e. additional 
interactions of the partons within the protons that un- 
derwent the hard interaction) are modelled using the 
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Herwig [11] and Jimmy [12] programs. For the con- 
struction of signal templates, the ti and single top 
quark production samples are generated for differ- 
ent assumptions on m top using six values (in GcV) 
namely (160,170,172.5,175,180,190), and with the 
largest samples at m top = 172.5 GeV. All it samples 
are normalised to the corresponding cross-sections, ob- 
tained with the latest theoretical computation approx- 
imating the NNLO prediction and implemented in the 
HATHOR package [13]. The predicted ti cross-section 
for a top quark mass of m t0 p = 172.5 GeV is 164.6 pb, 
with an uncertainty of about 8%. 

The production of W bosons or Z bosons in associ- 
ation with jets is simulated using the Alpgen genera- 
tor [14] interfaced to the Herwig and Jimmy packages. 
Diboson production processes (WW, WZ and ZZ) are 
produced using the Herwig generator. All Monte Carlo 
samples are generated with additional multiple soft 
proton-proton interactions. These simulated events are 
re-weighted such that the distribution of the number 
of interactions per bunch crossing (pileup) in the simu- 
lated samples matches that in the data. The mean num- 
ber of primary vertices per bunch crossing for the data 
of this analysis is about four. The samples are then pro- 
cessed through the GEANT4 [15] simulation [16] and 
the reconstruction software of the ATLAS detector. 



4 Event selection 

In the signal events the main reconstructed objects in 
the detector are electron and muon candidates as well 
as jets and missing transverse momentum (_E™ 1SS ). An 
electron candidate is defined as an energy deposit in the 
electromagnetic calorimeter with an associated well- 
reconstructed track. Electron candidates are required 
to have transverse energy Et > 25 GeV and | fluster | < 
2.47, where ry c i ustcr is the pseudorapidity of the electro- 
magnetic cluster associated with the electron. Candi- 
dates in the transition region between the barrel and 
end-cap calorimeter, i.e. candidates fulfilling 1.37 < 
hciustcrl < 1-52, are excluded. Muon candidates are 
reconstructed from track segments in different layers 
of the muon chambers. These segments are combined 
starting from the outermost layer, with a procedure 
that takes material effects into account, and matched 
with tracks found in the inner detector. The final candi- 
dates are refitted using the complete track information, 
and are required to satisfy pr > 20 GeV and \r]\ < 2.5. 
Isolation criteria, which restrict the amount of energy 
deposits near the candidates, are applied to both elec- 
tron and muon candidates to reduce the background 
from hadrons mimicking lepton signatures and back- 
grounds from heavy flavour decays inside jets. For elec- 



trons, the energy not associated to t he electron c luster 
and contained in a cone of AR — \j A<p + Aif = 0.2 
must not exceed 3.5 GeV, after correcting for energy de- 
posits from pileup, which in the order of 0.5 GeV. For 
muons, the sum of track transverse momenta and the 
total energy deposited in a cone of AR = 0.3 around 
the muon are both required to be less than 4 GeV. 

Jets are reconstructed with the anti-fc t algorithm [17] 
with R = 0.4, starting from energy clusters of adjacent 
calorimeter cells called topological clusters [18]. These 
jets are calibrated first by correcting the jet energy 
using the scale established for electromagnetic objects 
(EM scale) and then performing a further correction to 
the hadronic energy scale using correction factors, that 
depend on energy and 77, obtained from simulation and 
validated with data [19]. Jet quality criteria [20] are ap- 
plied to identify and reject jets reconstructed from en- 
ergies not associated to energy deposits in the calorime- 
ters originating from particles emerging from the bunch 
crossing under study. The jets failing the quality crite- 
ria, which may have been reconstructed from various 
sources such as calorimeter noise, non-collision beam- 
related background, and cosmic-ray induced showers, 
can efficiently be identified [20] . 

The reconstruction of E!f las is based upon the vec- 
tor sum of calorimeter energy deposits projected onto 
the transverse plane. It is reconstructed from topolog- 
ical clusters, calibrated at the EM scale and corrected 
according to the energy scale of the associated physics 
object. Contributions from muons are included by us- 
ing their momentum measured from the track and muon 
spectrometer systems in the E™ lss reconstruction. 

Muons reconstructed within a AR = 0.4 cone of 
a jet satisfying p T > 20 GeV are removed to reduce 
the contamination caused by muons from hadron de- 
cays within jets. Subsequently, jets within AR = 0.2 
of an electron candidate are removed to avoid double 
counting, which can occur because electron clusters are 
usually also reconstructed as jets. 

Reconstruction of top quark pair events is facilitated 
by the ability to tag jets originating from the hadroni- 
sation of 6-quarks. For this purpose, a neural-net-based 
algorithm [21], relying on vertex properties such as the 
decay length significance, is applied. The chosen work- 
ing point of the algorithm corresponds to a 6-tagging 
efficiency of 70% for jets originating from 6-quarks in 
simulated ti events and a light quark jet rejection factor 
of about 100. Irrespective of their origin, jets tagged by 
this algorithm are called 6-jets in the following, whereas 
those not tagged are called light jets. 

The signal is characterised by an isolated lepton 
with relatively high p^, E™ ss arising from the neutrino 
from the leptonic W boson decay, two 6-quark jets, and 
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Process 


Id-analysis 


2d-analysis 


e+jets 


//+jets 


e+jets 


/x+jets 


tt signal 

Single top (signal) 

Z+jets 

zz/wz/ww 

W+jets (data) 
QCD multijet (data) 


990 ± 40 
43 ± 2 


1450 ± 50 
53 ± 3 


3400 ± 200 


5100 + 300 


190 ± 10 

83 ± 8 

11± 2 

700 ± 500 

200 ± 200 


280+ 20 

100+ 8 

18+ 2 

1100 + 800 

400 ± 400 


12 ± 3 

2±<1 

80 ± 60 

50 ± 50 


8± 3 

2±<1 

100 ± 70 

40 ± 40 


Signal + background 


1180 ± 80 


1650 ± 80 


4500 ± 500 


6900 ± 900 


Data 


1151 


1724 


4556 


7225 



Table 1 The observed numbers of events in the data in the e+jets and /x+jets channels, for the two analyses after the common 
event selection and additional analysis-specific requirements. In addition, the expected numbers of signal and background events 
corresponding to the integrated luminosity of the data are given, where the single top quark production events are treated as 
signal for the Id-analysis, and as background for the 2d-analysis. The Monte Carlo estimates assume SM cross-sections. The 
Vy+jets and QCD multijet background contributions are estimated from ATLAS data. The uncertainties for the estimates 
include different components detailed in the text. All predicted event numbers are quoted using one significant digit for the 
uncertainties, i.e. the trailing zeros are insignificant. 



two light quark jets from the hadronic W boson decay. 
The selection of events consists of a series of require- 
ments on general event quality and the reconstructed 
objects designed to select the event topology described 
above. The following event selections are applied: 

— it is required that the appropriate single electron 
or single muon trigger has fired (with thresholds at 
20 GeV and 18 GeV, respectively); 

— the event must contain one and only one recon- 
structed lepton with £+ > 25 GeV for electrons 
and pt > 20 GeV for muons which, for the e+jets 
channel, should also match the corresponding trig- 
ger object; 

— in the /z+jets channel, E™ lss > 20 GeV and in addi- 
tion E™ lss + m^ > 60 GeV is required 2 ; 

— in the e+jets channel more stringent cuts on E!f lsa 



and 



'w 



are required because of the higher level 



of QCD multijet background, these being E, 



miss 
T 



> 



35 GeV and 



'w 



> 25 GeV; 



— the event is required to have > 4 jets with p^ > 
25 GeV and |?7| < 2.5. It is required that at least 
one of these jets is a 6-jet. 

This common event selection is augmented by ad- 
ditional analysis-specific event requirements described 
next. 



5 Specific event requirements 

To optimise the expected total uncertainty on m t0 p, 
some specific requirements are used in addition to the 
common event selection. 

2 Here rn^ is the VK-boson transverse mass, defined as 
•^/2pT,f Pt,v [1 — cos(<^ — 4> v )], where the measured E!^- laa 
vector provides the neutrino (u) information. 



For the Id-analysis, three additional requirements 
are applied. Firstly, only events with a converging like- 
lihood fit (see Section 6) with a logarithm of the likeli- 
hood value In L > —50 are retained. Secondly, all jets in 
the jet triplet assigned to the hadronic decay of the top 
quark are required to fulfill pt > 40 GeV, and thirdly 
the reconstructed W boson mass must lie within the 
range 60 GeV - 100 GeV. 

For the 2d-analysis the additional requirement is 
that only light jet pairs (see Section 7) with an invari- 
ant mass in the range 50 GeV -110 GeV are considered 
for the x 2 nt - 

The numbers of events observed and expected, with 
the above selection and these additional analysis-specific 
requirements, are given in Table 1 for both channels and 
both analyses. For all Monte Carlo estimates, the un- 
certainties are the quadratic sum of the statistical un- 
certainty, the uncertainty on the 6-tagging efficiencies, 
and a 3.7% uncertainty on the luminosity [22,23]. For 
the QCD multijet and the VF+jets backgrounds, the 
systematic uncertainty estimated from data [24] domi- 
nates and is used instead. 

For both analyses and channels, the observed distri- 
butions for the leptons, jets, and kinematic properties 
of the top quark candidates such as their transverse 
momenta, are all well-described by the sum of the sig- 
nal and background estimates. This is demonstrated for 
the properties of the selected jets, before applying the 
analysis specific requirements, for both channels in Fig- 
ure 1. The jet multiplicities, shown in Figure l(a, b), 
as well as the distributions of kinematic properties of 
jets like transverse momenta, Figure l(c, d), and the 
77 distributions, Figure l(e, f), are all well-described 
within the uncertainty band of the prediction. The size 
of the uncertainty band is dominated by the uncertain- 
ties on the background contributions estimated from 
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(e) e+jets channel 
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Fig. 1 Distributions for the selected events of the common event selection in the e+jets channel on the left and the /x+jets 
channel on the right. Shown are (a, b) the measured jet multiplicities, (c, d) the px, and (e, f) the r\ distributions of all selected 
jets. The hatched area is the total uncertainty on the prediction described in the text. In (c, d) the rightmost bin also contains 
the overflow. 
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data. The largest differences between the central values 
of the combined prediction and the data is observed for 
the rapidity distribution, with the data being higher, 
especially at central rapidities. Based on the selected 
events, the top quark mass is measured in two ways as 
described below. 



6 The Id-analysis 

The Id-analysis is a one-dimensional template analysis 
using the reconstructed mass ratio: 



R 



32 



"top 



m 



w 



Here m\^° and m^f" are the per event reconstructed 
invariant masses of the hadronically decaying top quark 
and W boson, respectively. 

To select the jet triplet for determining the two 
masses, this analysis utilises a kinematic fit maximis- 
ing an event likelihood. This likelihood relates the ob- 
served objects to the tt decay products (quarks and 
leptons) predicted by the NLO signal Monte Carlo, al- 
beit in a Leading Order (LO) kinematic approach, using 
ti — > ivbi gi^^had- In this procedure, the measured jets 
relate to the quark decay products of the W boson, gi 
and (72, and to the 6-quarks, bi and frhad, produced in 
the top quark decays. The _E™ 1SS vector is identified with 
the transverse momentum components of the neutrino, 
p x<v and p y<v . 

The likelihood is defined as a product of transfer 
functions (T), Breit-Wigner (B) distributions, and a 
weight Wbtag accounting for the 6-tagging information: 



L = T[E ]ctl 



r{E jct2 

TiEl 



E h 



T E ]ct \E, 



T (E™ 1SS 



T[E icti \E q2 yT{Er s \p x 
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m{qiq 2 b had )\m t ^ ,Ft 



'top 



top 



B 



/ fj t m rcco.likc 7-1 

m(tvbz)\m top ' ,r top 



w ht 



The generator predicted quantities are marked with 
a circumflex (e.g. Eh had ), i.e. the energy of the &-quark 
from the hadronic decay of the top quark. The quanti- 
ties my/ and Fw (which amounts to about one fifth of 
the Gaussian resolution of the m^° distribution) are 
taken from Ref. [6], and m™p°' ' ° is the likelihood esti- 
mator for the top quark mass, i.e. the per event result 
of maximising this likelihood. Transfer functions are de- 
rived from the MC@NLO tt signal Monte Carlo sample 
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Fig. 2 Id-analysis: Performance of the likelihood fit in the 
e+jets channel. Shown in (a) are the predicted In L distri- 
butions for various jet permutations in the tt signal Monte 
Carlo. The figures (b, c) compare two output variables of the 
likelihood fit as observed in the data with their respective 
prediction. These are (b) the In L value, and (c) the px of 
the b-jet associated to the hadronic decay of the top quark. 



Measurement of the Top Quark Mass with the Template Method using ATLAS Data 



at an input mass of mtop = 172.5 GeV, based on recon- 
structed objects that are matched to their generator 
predicted quarks and leptons. When using a maximum 
separation of AR = 0.4 between a quark and the corre- 
sponding jet, the fraction of events with four matched 
jets from all selected events amounts to 30% - 40%. The 
transfer functions are obtained in three bins of r\ for the 
energies of 6-quark jets, £j e t, and -Ejet, , light quark jets, 
£j e t and E^ , the energy, E e , (or transverse momen- 
tum, pr,fj.) of the charged lepton, and the two compo- 
nents of the £ T liss , Ef lss and E™ [ss . In addition, the 
likelihood exploits the values of my/ and TV 1° con- 
strain the reconstructed leptonic, m(£v), and hadronic, 
m(qi (72)7 W boson masses using Breit-Wigner distribu- 
tions. Similarly, the reconstructed leptonic, m(£vbi), 
and hadronic, m{q\ q% &had), top quark masses are con- 
strained to be identical, where the width of the cor- 
responding Breit Wigner distribution is identified with 
the predicted r t0 p (using its top quark mass depen- 
dence) [6] . Including the ^-tagging information into the 
likelihood as a weight Wbtagj derived from the efficiency 
and mistag rate of the 6-tagging algorithm, and as- 
signed per jet permutation according to the role of each 
jet for a given jet permutation, improves the selection 
of the correct jet permutation. As an example, for a 
permutation with two 6-jets assigned to the 6-quark po- 
sitions and two light jets to the light quark positions, 
the weight Wbtag amounts to 0.48, i.e. it corresponds to 
the square of the 6-tagging efficiency times the square 
of one minus the fake rate, both given in Section 4. 

With this procedure, the correct jet triplet for the 
hadronic top quark is chosen in about 70% of simu- 
lated signal events with four matched jets. However, if 
i? 32 from the likelihood fit, i.e. calculated from m^™' 
and m\^°' c , is taken, a large residual jet energy scale 
(JES) dependence of R32 remains. This is because in 
the fit m^ co is constrained to m\y, while m{£ p ° is only 
constrained to be equal for the leptonic and hadronic 
decays of the top quarks. This spoils the desired event- 
by-event reduction of the JES uncertainty in the ratio 
-R32 [25]. To make best use of the high selection effi- 
ciency for the correct jet permutation from the likeli- 
hood fit, and the stabilisation of R32 against JES varia- 
tions, the jet permutation derived in the fit is used, but 



m\\'". m;.™ and therefore R32, are constructed from 



,rcco 
l W > 



the unconstrained four- vectors of the jet triplet as given 
by the jet reconstruction. 

The performance of the algorithm, shown in Fig- 
ure 2 for the e+jets channel, is similar for both chan- 
nels. The likelihood values of wrong jet permutations 
for signal events from the large MC@NLO sample are 
frequently considerably lower than the ones for the cor- 
rect jet permutations, as seen in Figure 2(a). For exam- 
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Fig. 3 Id-analysis: The reconstructed R32 constructed 
from the selected jet permutation using the unconstrained 
four-vectors of the jet triplet for (a) the e+jets channel, and 
(b) the /x+jets channel. The rightmost bins also contain the 
overflow. 



pie, the distribution for the jet permutation in which 
the jet from the 6-quark from the leptonically decaying 
top quark is exchanged with one light quark jet from 
the hadronic W boson decay has a second peak at about 
ten units lower than the one for the correct jet permuta- 
tion. The actual distribution of In L values observed in 
the data is well-described by the signal plus background 
predictions, as seen in Figure 2(b). The kinematic dis- 
tributions of the variables used in the transfer functions 
are also well-described by the predictions, as shown in 
Figure 2(c), for the example of the resulting px of the 
6-jet associated to the hadronic decay of the top quark. 
The resulting R32 distributions for both channels are 
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shown in Figure 3. They are also well accounted for by 
the predictions. 

Signal templates are derived for the -R32 distribu- 
tion for all m t0 p dependent samples, consisting of the 
ti signal events, together with single top quark produc- 
tion events. This procedure is adopted, firstly, because 
single top quark production, although formally a back- 
ground process, still carries information about the top 
quark mass and, secondly, by doing so m to p indepen- 
dent background templates can be used. The templates 
are constructed for the six m top choices using the specif- 
ically generated Monte Carlo samples, see Section 3. 

The R32 templates are parameterised with a func- 
tional form given by the sum of a ratio of two corre- 
lated Gaussians and a Landau function. The ratio of 
two Gaussians [26] is motivated as a representation of 
the ratio of two correlated measured masses. The Lan- 
dau function is used to describe the tails of the dis- 
tribution stemming mainly from wrong jet-triplet as- 
signments. The correlation between the two Gaussian 
distributions is fixed to 50%. A simultaneous fit to all 
templates per decay channel is used to derive a continu- 
ous function of m t0 p that interpolates the R32 shape dif- 
ferences among all mass points with mtop in the range 
described above. This approach rests on the assump- 
tion that each parameter has a linear dependence on 
the top quark mass, which has been verified for both 
channels. The fit minimises a \ 2 built from the -R32 
distributions at all mass points simultaneously. The x 2 
is the sum over all bins of the difference squared be- 
tween the template and the functional form, divided by 
the statistical uncertainty squared in the template. The 
combined fit adequately describes the R32 distributions 
for both channels. In Figure 4(a) the sensitivity to m top 
is shown in the e+jets channel by the superposition of 
the signal templates and their fits for four of the six 
input top quark masses assumed in the simulation. 

For the background template, the mtop independent 
parts, see Table 1, are treated together. Their individ- 
ual distributions, taken either from Monte Carlo or data 
estimates as detailed above, are summed, and a Lan- 
dau distribution is chosen to parameterise their R32 
distribution. For each channel this function adequately 
describes the background distribution as shown in Fig- 
ure 4(b) for the e+jets channel, which has a larger back- 
ground contribution than the /^+jets channel. 

Signal and background probability density functions, 
P s ig(R32\mto P ) and -P D k g (-R32), respectively, are used in 
a binned likelihood fit to the data using a number of 
bins, -/Vbins- The likelihood reads: 

£(i?32|m to p) = £sha,pe(R32\mto P ) X £bkg(-R32) , 
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Fig. 4 Id-analysis: Template parameterisations for (a) sig- 
nal and (b) background contributions in the e+jets channel. 
The background fit is labelled Pbkg- 
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with: 



Ai = (N - n hkg ) ■ P s i g (i?32|m t op) l + n hkg ■ P hkg (R 3 2)i 



N h 



N = Y, N i = ™sig + 



"bkg 



The variable N[ denotes the number of events observed 
per bin, and n s j g and ribkg denote the total numbers 
of signal and background events to be determined. The 
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term £ s hape accounts for the shape of the R32 distribu- 
tion and its dependence on the top quark mass mtop- 
The term £bkg constrains the total number of back- 
ground events, nbkg, using its prediction, n^™ , and 
the background uncertainty, chosen to be 50%, see Ta- 
ble 1. In addition, the number of background events 
is restricted to be positive. The two free parameters 
of the fit are the total number of background events, 
nbkg, and mtop- The performance of this algorithm is as- 
sessed with the pseudo-experiment technique. For each 
m-top value, distributions from pseudo-experiments are 
constructed by random sampling of the simulated sig- 
nal and background events used to construct the corre- 
sponding templates. Using Poisson statistics, the num- 
bers of signal events and total background events in 
each pseudo-experiment are fluctuated around the ex- 
pectation values, either calculated assuming SM cross- 
sections and the integrated luminosity of the data, or 
taken from the data estimate. A good linearity is found 
between the input top quark mass used to perform the 
pseudo-experiments, and the result of the fit. Within 
their statistical uncertainties, the mean values and width 
of the pull distributions are consistent with the ex- 
pectations of zero and one, respectively. The expected 
statistical uncertainties (mean ± RMS) obtained from 
pseudo-experiments with an input top quark mass of 
m-top = 172.5 GeV, and for a luminosity of 1 fb -1 , are 
1.36±0.16 GeV and 1.11 ±0.06 GeV for the e+jets and 
/x+jets channels, respectively. 



7 The 2d-analysis 

In the 2d-analysis, similarly to Ref. [27], m t0 p and a 
global jet energy scale factor (JSF) are determined 
simultaneously by using the m\^° and m^° distri- 
butions 3 . Instead of stabilising the estimator of m-top 
against JES variations as done for the Id-analysis, the 
emphasis here is on an in-situ jet scaling. A global JSF 
(averaged over 77 and pr) is obtained, which is mainly 
based on the observed differences between the predicted 
m!^ co distribution and the one observed for the data. 
This algorithm predicts which global JSF correction 
should be applied to all jets to best fit the data. Due 
to this procedure, the JSF is sensitive not only to the 
JES, but also to all possible differences in data and 
predictions from specific assumptions made in the sim- 
ulation that can lead to differences in the observed jets. 
These comprise: the fragmentation model, initial state 
and final state QCD radiation (ISR and FSR), the un- 

3 Although for the two analyses m^p and rn^ c ° are calcu- 
lated differently, the same symbols are used to indicate that 
these are estimates of the same quantities. 



derlying event, and also pileup. In this method, the sys- 
tematic uncertainty on mt op stemming from the JES is 
reduced and partly transformed into an additional sta- 
tistical uncertainty on m top due to the two-dimensional 
fit. The precisely measured values of mw and TV [6] 
are used to improve on the experimental resolution of 
ml° c ° by relating the observed jet energies to the cor- 
responding parton energies as predicted by the signal 
Monte Carlo (i.e. to the two quarks from the hadronic 
W boson decay, again using LO kinematics). Thereby, 
this method offers a complementary determination of 
m top to the Id-analysis method, described in Section 6, 
with different sensitivity to systematic effects and data 
statistics. 

For the events fulfilling the common requirements 
listed in Section 4, the jet triplet assigned to the 
hadronic top quark decay is constructed from any &-jet, 
together with any light jet pair with a reconstructed 
m^ co within 50 GeV - 110 GeV. Amongst those, the 
jet triplet with maximum px is chosen as the top quark 
candidate. For the light jet pair, i.e. for the hadronic 
W boson decay candidates, a kinematic fit is then per- 
formed by minimising the following \ 2 : 



X 



E 



-Ejet,i(l 



0"(-Ejet,i) 



-Mjet,jet(Qi,a 2 ) -m w 



with respect to parton scale factors (ai) for the jet en- 
ergies. The \ 2 comprises two components. The first 
component is the sum of squares of the differences 
of the measured and fitted energies of the two recon- 
structed light jets, -Ejct.i, individually divided by the 
squares of their p^- and 77-dependent resolutions ob- 
tained from Monte Carlo simulation, cr(-Ej eti i). The sec- 
ond term is the difference of their two-jet invariant 
mass, Mjotjot, and mw, divided by the W boson width. 
From these jets the two observables m^ co and m^p 



are constructed. The 



'w 



is calculated using the re- 



constructed light jet four-vectors (i.e. jet energies are 
not corrected using a\), retaining the full sensitivity of 



'w 



to the JSF. In contrast, 



Hop 



is calculated from 



these light jet four-vectors scaled to the parton level 
(i.e. jet energies are corrected using ai) and the above 
determined 6-jet. In this way light jets in m^p exhibit a 
much reduced JES sensitivity by construction, and only 
the 6-jet is directly sensitive to the JES. The m^ co and 
m\^° distributions are shown in Figure 5 for both lep- 
ton channels, together with the predictions for signal 
and background. These, in both cases describe the ob- 
served distributions well. The correlation of these two 
observables is found to be small for data and predic- 
tions, and amounts to about —0.06. 

Templates are constructed for m\ c ^° as a function 
of an input top quark mass in the range 160 GeV - 
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and mj™°, observed in the data together with the 
signal and background predictions. Shown are (a, c) the e+jets channel, and (b, d) the /x+jets channel. 



190 GeV, and of an input value for the JSF in the 
range 0.9 - 1.1, and, finally, for m^p° as a function 
of the assumed JSF for the same range. The signal 
templates for the to^ co and m}^ distributions, shown 
for the /x+jets channel and for JSF=1 in Figure 6(a) 
and 6(b), are fitted to a sum of two Gaussian func- 
tions for m^ co , and to the sum of a Gaussian and a 



Landau function for 



"top 



Since, for this analysis, the 



background templates are constructed including single 
top quark production events, the background fit for the 



m\^° distribution is assumed to be mt op dependent. For 
the background, the to^ co distribution, again shown for 
the /i+jets channel in Figure 6(c), is fitted to a Gaus- 
sian function and the m^" distribution, Figure 6(d), to 
a Landau function. For all parameters of the functions 
that also depend on the JSF, a linear parameterisation 
is chosen. The quality of all fits is good for the signal 
and background contributions and for both channels. 

Signal and background probability density functions 
for the m^p° and m^° distributions are used in an 
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Fig. 6 2d-analysis: Template parameterisations for the m^ co and m rec 



,.,,, distributions for signal and background events 
for the /i+jets channel. Shown are (a, b) the m^ co and mj^p distributions for signal events, and (c, d) the corresponding 
distributions for background events, see Table 1. All distributions are for JSF=1. 



unbinncd likelihood fit to the data for all events, i — 
1, . . . N. The likelihood function maximised is: 

Ahapc(m^ co , m™ p |m t0 p, JSF, n bkg ) = 

N 

Y[Ptop(ml™\mto P , JSF, n b k g ) % x 
»=i 

iV(n*w°| JSF, nbk g )» , 

with: 

Pto P = {N- n bkg ) ■ P£ s p (mZp°\mtop, JSF), + 

n bkg .P^ k p g K-°|m top ,JSF),, 
P W = (N- n bkg ) ■ P^ s (m- co | JSF), + 

n bkg • P^WIJSF)* . 

The three parameters to be determined by the fit are 
m t0 p, the JSF and n bkg . Using pseudo-experiments, a 
good linearity is found between the input top quark 
mass used to perform the pseudo-experiments, and the 
result of the fits. The residual dependence of the recon- 
structed m top is about 0.1 GeV for a JSF shift of 0.01 



for both channels, which results in a residual system- 
atic uncertainty due to the JES. Within their statis- 
tical uncertainties, the mean values and widths of the 
pull distributions are consistent with the expectations 
of zero and one, respectively. Finally, the expected sta- 
tistical plus JSF uncertainties (mean ± RMS) obtained 
from pseudo-experiments at an input top quark mass 
of mtop = 172.5 GeV, and for a luminosity of 1 fb -1 , 
are 1.20 ±0.08 GeV and 0.94 ±0.04 GeV for the e+jets 
and /z+jets channel, respectively. 



8 Top quark mass measurement 

8.1 Evaluation of systematic uncertainties 

Each source of uncertainty considered is investigated, 
when possible, by varying the respective quantities 
by ±1<T with respect to the default value. Using 
the changed parameters, pseudo-experiments are either 
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performed directly or templates are constructed and 
then used to generate pseudo-experiments, without al- 
tering the probability density function parameterisa- 
tions. The difference of the results for m top compared 
to the standard analysis is used to determine the sys- 
tematic uncertainties. For the 2d-analysis, in any of the 
evaluations of the systematic uncertainties, apart from 
the JES variations, the maximum deviation of the JSF 
from its nominal fitted value is ±2.5%. 

All sources of systematic uncertainties investigated, 
together with the resulting uncertainties, are listed in 
Table 2. The statistical precision on m top obtained 
from the Monte Carlo samples is between 0.2 GeV and 
0.5 GeV, depending on the available Monte Carlo statis- 
tics. For some sources, pairs of statistically independent 
samples are used. For other sources, the same sample 
is used, but with a changed parameter. In this case the 
observed mtop values for the central and the changed 
sample arc statistically highly correlated. In all cases, 
the actual observed difference is quoted as the system- 
atic uncertainty on the corresponding source, even if it 
is smaller than the statistical precision of the difference. 
The total uncertainty is calculated as the quadratic sum 
of all individual contributions, i.e. neglecting possible 
correlations. The estimation of the uncertainties from 
the individual contributions is described in the follow- 
ing. 

Jet energy scale factor: This is needed to sepa- 
rate the quoted statistical uncertainty on the result of 
the 2d-analysis into a purely statistical component on 
TO top analogous to the one obtained in an Id-analysis, 
and the contribution stemming from the simultaneous 
determination of the JSF. This uncertainty is evalu- 
ated for the 2d-analysis by in addition performing a 
one-dimensional (i.e. JSF-constraint) fit to the data, 
with the JSF fixed to the value obtained in the two- 
dimensional fit. The quoted statistical precision on mtop 
is the one from the one-dimensional fit. The contribu- 
tion of the JSF is obtained by quadratically subtract- 
ing the statistical uncertainties on m t0 p for the one- 
dimensional and two-dimensional fit of the 2d-analysis. 

Method calibration: The limited statistics of the 
Monte Carlo samples leads to a systematic uncertainty 
in the template fits, which is reflected in the residual 
mass differences between the fitted and the input mass 
for a given Monte Carlo sample. The average difference 
observed in the six samples with different input masses 
is taken as the uncertainty from this source. 

Signal Monte Carlo generator: The systematic 
uncertainty related to the choice of the generator pro- 
gram is accounted for by comparing the results of pseudo- 
experiments performed with either the MC@NLO or 



the Powheg samples [28] both generated with mtop = 
172.5 GeV. 

Hadronisation: Signal samples for m t0 p = 
172.5 GeV from the Powheg event generator are pro- 
duced with either the Pythia [29] or Herwig [11] 
program performing the hadronisation. One pseudo- 
experiment per sample is performed and the full dif- 
ference of the two results is quoted as the systematic 
uncertainty. 

Pileup: To investigate the uncertainty due to ad- 
ditional proton-proton interactions which may affect 
the jet energy measurement, on top of the component 
that is already included in the JES uncertainty dis- 
cussed below, the fit is repeated in data and simulation 
as a function of the number of reconstructed vertices. 
Within statistics, the measured mtop is independent of 
the number of reconstructed vertices. This is also ob- 
served when the data are instead divided into data peri- 
ods according to the average numbers of reconstructed 
vertices. In this case, the subsets have varying contri- 
butions from pileup from preceding events. 

However, the effect on m t0 p due to any residual 
small difference between data and simulation in the 
number of reconstructed vertices was assessed by com- 
puting the weighted sum of a linear interpolation of 
the fitted masses as a function of the number of pri- 
mary vertices. In this sum the weights are the relative 
frequency of observing a given number of vertices in the 
respective sample. The difference of the sums in data 
and simulation is taken as the uncertainty from this 
source. 

Underlying event: This systematic uncertainty is 
obtained by comparing the AcerMC [30, 31] central 
value, defined as the average of the highest and the low- 
est masses measured on the ISR/FSR variation samples 
described below, with a dataset with a modified under- 
lying event. 

Colour reconnection: The systematic uncertainty 
due to colour reconnection is determined using Ac- 
erMC with Pythia with two different simulations of 
the colour reconnection effects as described in Refs. [32- 
34]. In each case, the difference in the fitted mass be- 
tween two assumptions on the size of colour reconnec- 
tion was measured. The maximum difference is taken as 
the systematic uncertainty due to colour reconnection. 

Initial and final state QCD radiation: Differ- 
ent amounts of initial and final state QCD radiation 
can alter the jet energies and the jet multiplicity of the 
events with the consequence of introducing distortions 
into the measured ml° p ° and m^, co distributions. This 
effect is evaluated by performing pseudo-experiments 
for which signal templates are derived from seven ded- 
icated AcerMC signal samples in which Pythia pa- 
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Id-analysis 


2d-analysis 


Combinations 


Correlation 




e+jets 


/i+jets 


e+jets 


^t+jets 


Id 


2d 


P 


Measured value of mtop 


172.93 


175.54 


174.30 


175.01 


174.35 


174.53 




Data statistics 


1.46 


1.13 


0.83 


0.74 


0.91 


0.61 




Jet energy scale factor 


na 


na 


0.59 


0.51 


na 


0.43 





Method calibration 


0.07 


< 0.05 


0.10 


< 0.05 


< 0.05 


0.07 





Signal MC generator 


0.81 


0.69 


0.39 


0.22 


0.74 


0.33 


1 


Hadronisation 


0.33 


0.52 


0.20 


0.06 


0.43 


0.15 


1 


Pileup 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


< 0.05 


<0.05 


1 


Underlying event 


0.06 


0.10 


0.42 


0.96 


0.08 


0.59 


1 


Colour reconnection 


0.47 


0.74 


0.32 


1.04 


0.62 


0.55 


1 


ISR and FSR (signal only) 


1.45 


1.40 


1.04 


0.95 


1.42 


1.01 


1 


Proton PDF 


0.22 


0.09 


0.10 


0.10 


0.15 


0.10 


1 


W+jets background normalisation 


0.16 


0.19 


0.34 


0.44 


0.18 


0.37 


1 


VF+jets background shape 


0.11 


0.18 


0.07 


0.22 


0.15 


0.12 


1 


QCD multijet background normalisation 


0.07 


< 0.05 


0.25 


0.33 


< 0.05 


0.20 


(1) 


QCD multijet background shape 


0.14 


0.12 


0.38 


0.30 


0.09 


0.27 


(1) 


Jet energy scale 


1.21 


1.25 


0.63 


0.71 


1.23 


0.66 


1 


6-jet energy scale 


1.09 


1.21 


1.61 


1.53 


1.16 


1.58 


1 


6-tagging efficiency and mistag rate 


0.21 


0.13 


0.31 


0.26 


0.17 


0.29 


1 


Jet energy resolution 


0.34 


0.38 


0.07 


0.07 


0.36 


0.07 


1 


Jet reconstruction efficiency 


0.08 


0.11 


< 0.05 


< 0.05 


0.10 


<0.05 


1 


Missing transverse momentum 


<0.05 


<0.05 


0.12 


0.16 


<0.05 


0.13 


1 


Total systematic uncertainty 


2.46 


2.56 


2.31 


2.57 


2.50 


2.31 




Total uncertainty 


2.86 


2.80 


2.46 


2.68 


2.66 


2.39 





Table 2 The measured values of m top and the contributions of various sources to the uncertainty of m top (in GeV) together 
with the assumed correlations p between analyses and lepton channels. Here '0' stands for uncorrelated, '1' for fully correlated 
between analyses and lepton channels, and '(1)' for fully correlated between analyses, but uncorrelated between lepton channels. 
The abbreviation 'na' stands for not applicable. The combined results described in Section 8.2 are also listed. 



rameters that control the showering are varied in ranges 
that are compatible with those used in the Perugia 
Hard/Soft tune variations [32]. The systematic uncer- 
tainty is taken as half the maximum difference between 
any two samples. 

Using different observables, the additional jet activ- 
ity accompanying the jets assigned to the top quark de- 
cays has been studied. For events in which one (both) 
W bosons from the top quark decays themselves de- 
cay into a charged lepton and a neutrino, the recon- 
structed jet multiplicities [35] (the fraction of events 
with no additional jet above a certain transverse mo- 
mentum [36]) are measured. The analysis of the recon- 
structed jet multiplicities is not sufficiently precise to 
constrain the presently used variations of Monte Carlo 
parameters. In contrast, for the ratio analysis [36] the 
spread of the predictions caused by the presently per- 
formed ISR variations is significantly wider than the 
uncertainty of the data, indicating that the present ISR 
variations are generous. 

Proton PDF: The signal samples are generated us- 
ing the CTEQ 6.6 [10] proton parton distribution func- 
tions, PDFs. These PDFs, obtained from experimental 
data, have an uncertainty that is reflected in 22 pairs 
of additional PDF sets provided by the CTEQ group. 
To evaluate the impact of the PDF uncertainty on the 
signal templates, the events are re-weighted with the 



corresponding ratio of PDFs, and 22 pairs of additional 
signal templates are constructed. Using these templates 
one pseudo-experiment per pair is performed. The un- 
certainty is calculated as half the quadratic sum of dif- 
ferences of the 22 pairs as suggested in Ref. [37]. 

W+jets background normalisation: The uncer- 
tainty on the IF+jets background determined from data 
is dominated by the uncertainty on the heavy flavour 
content of these events and amounts to ±70%. The dif- 
ference in mtop obtained by varying the normalisation 
by this amount is taken as the systematic uncertainty. 

W+jets background shape: The impact of the 
variation of the shape of the 1^+jets background con- 
tribution is studied using a re-weighting algorithm [24] 
which is based on changes observed on stable parti- 
cle jets when model parameters in the Alpgen Monte 
Carlo program are varied. 

QCD multijet background normalisation: The 

estimate for the background from QCD multijet events 
determined from data is varied by ±100% to account for 
the current understanding of this background source [24] 
for the signal event topology. 

QCD multijet background shape: The uncer- 
tainty due to the QCD background shape has been es- 
timated comparing the results from two data driven 
methods, for both channels, see Ref. [24] for details. 
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For this uncertainty pseudo-experiments are performed 
on QCD background samples with varied shapes. 

Jet energy scale: The jet energy scale is derived 
using information from test-beam data, LHC collision 
data and simulation. Since the energy correction pro- 
cedure involves a number of steps, the JES uncertainty 
has various components originating from the calibration 
method, the calorimeter response, the detector sim- 
ulation, and the specific choice of parameters in the 
physics model employed in the Monte Carlo event gen- 
erator. The JES uncertainty varies between ±2.5% and 
±8% in the central region, depending on jet px and 
7/ as given in Ref. [19]. These values include uncer- 
tainties in the flavour composition of the sample and 
mis-measurements from jets close by. Pileup gives an 
additional uncertainty of up to ±2.5% (±5%) in the 
central (forward) region. Due to the use of the observ- 
able -R32 for the Id-analysis, and to the simultaneous fit 
of the JSF and mt op for the 2d-analysis, which mitigate 
the impact of the JES on TO top differently, the system- 
atic uncertainty on the determined mtop resulting from 
the uncertainty of the jet energy scale is less than 1%, 
i.e. much smaller than the JES uncertainty itself. 

Relative b-jet energy scale: This uncertainty is 
uncorrelated with the jet energy scale uncertainty and 
accounts for the remaining differences between jets orig- 
inating from light quarks and those from 6-quarks af- 
ter the global JES has been determined. For this, an 
extra uncertainty ranging from ±0.8% to ±2.5% and 
depending on jet pt and rj is assigned to jets arising 
from the fragmentation of 6-quarks, due to differences 
between light jets and gluon jets, and jets containing 
6-hadrons. This uncertainty decreases with px, and the 
average uncertainty for the spectrum of jets selected in 
the analyses is below ±2%. 

This additional systematic uncertainty has been ob- 
tained from Monte Carlo simulation and was also veri- 
fied using &-jets in data. The validation of the 6-jet en- 
ergy scale uncertainty is based on the comparison of the 
jet transverse momentum as measured in the calorime- 
ter to the total transverse momentum of charged par- 
ticle tracks associated to the jet. These transverse mo- 
menta are evaluated in the data and in Monte Carlo 
simulated events for inclusive jet samples and for 6-jet 
samples [19]. Moreover, the jet calorimeter response un- 
certainty has been evaluated from the single hadron re- 
sponse. Effects stemming from 6-quark fragmentation, 
hadronisation and underlying soft radiation have been 
studied using different Monte Carlo event generation 
models [19]. 

b-tagging efficiency and mistag rate: The 
^-tagging efficiency and mistag rates in data and Monte 
Carlo simulation are not identical. To accommodate 



this, 6-tagging scale factors, together with their uncer- 
tainties, are derived per jet [21,38]. They depend on the 
jet px a< n d f] and the underlying quark-flavour. For the 
default result the central values of the scale factors are 
applied, and the systematic uncertainty is assessed by 
changing their values within their uncertainties. 

Jet energy resolution: To assess the impact of 
this uncertainty, before performing the event selection, 
the energy of each reconstructed jet in the simulation is 
additionally smeared by a Gaussian function such that 
the width of the resulting Gaussian distribution corre- 
sponds to the one including the uncertainty on the jet 
energy resolution. The fit is performed using smeared 
jets and the difference to the default m tov measurement 
is assigned as a systematic uncertainty. 

Jet reconstruction efficiency: The jet reconstruc- 
tion efficiency for data and the Monte Carlo simulation 
are found to be in agreement with an accuracy of better 
than ±2% [19]. To account for this, jets are randomly 
removed from the events using that fraction. The event 
selection and the fit are repeated on the changed sam- 
ple. 

Missing transverse momentum: The E™ ss is 
used in the event selection and also in the likelihood for 
the Id-analysis, but is not used in the 77i top estimator 
for either analysis. Consequently, the uncertainty due to 
any mis-calibration is expected to be small. The impact 
of a possible mis-calibration is assessed by changing the 
measured E™ lss within its uncertainty. 

The resulting sizes of all uncertainties are given in 
Table 2. They are also used in the combination of results 
described below. The three most important sources of 
systematic uncertainty for both analyses are the rel- 
ative 6-jet to light jet energy scale, the modelling of 
initial and final state QCD radiation, and the light jet 
energy scale. Their impact on the precision on m top are 
different as expected from the difference in the estima- 
tors used by the two analyses. 



8.2 Results 

Figure 7 shows the results of the Id-analysis when per- 
formed on data. For both channels, the fit function de- 
scribes the data well, with a x 2 / do1 of 21 / 23 (39/23) 
for the e+jcts (/j+jets) channels. The observed statis- 
tical uncertainties in the data are consistent with the 
expectations given in Section 6 with the e±jets chan- 
nel uncertainty being slightly higher than the expected 
uncertainty of 1.36 ± 0.16 GeV. The results from both 
channels are statistically consistent and are: 

m top = 172.9 ± 1.5 sta t ± 2.5 syst GeV (Id e+jets), 
m t0 p = 175.5 ± l.l s tat ± 2.6 syst GeV (Id /i±jets). 
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Fig. 7 Id-analysis: The R32 distribution observed in the 
data together with the signal and background contributions 
determined by the fit. The distributions are for (a) the e+jets 
channel and (b) the /x+jets channel. The data points are 
shown with their statistical uncertainties. 

Figure 8 shows the results of the 2d-analysis when 
performed on data for the e+jets and /x+jets channels. 
Again the fit functions describe the observed distribu- 
tions well, with a x 2 / dor of 47/38 (51/38) for the sum of 
the m^- co and m^™ distributions in the e+jets (/x+jets) 
channels. The two-dimensional uncertainty ellipses for 
both channels are shown in Figure 9. The results from 
both channels are: 

m top = 174.3 ± I.0 s tat+JSF ± 2.2 syst GeV (2d e+jets), 
m top = 175.0 ± 0.9 stat+ jsF ± 2.5 syst GeV (2d xt+jcts). 

Within statistical uncertainties these results are consis- 
tent with each other, and the observed statistical un- 
certainties in the data are in accord with the expec- 
tations given in Section 7, however, for this analysis, 
with the e+jets channel uncertainty being slightly lower 
than the expected uncertainty of 1.20 ± 0.08 GeV. The 
corresponding values for the JSF are 0.985 ± 0.008 and 
0.986 ±0.006 in the e+jets and /x+jets channels, respec- 



tively, where the uncertainties are statistical only. The 
JSF values fitted for the two channels are consistent 
within their statistical uncertainty. For both channels, 
the correlation of m t0 p and the JSF in the fits is about 
-0.57. 

When separating the statistical and JSF component 
of the result as explained in the discussion of the JSF 
uncertainty evaluation in Section 8.1, the result from 
the 2d-analysis yields: 

m top = 174.3 ± 0.8 sta t ± 2.3 syst GeV (2d e+jets), 
m top = 175.0 ± 0.7 stat ± 2.6 syst GeV (2d /x+jcts). 

These values together with the breakdown of uncertain- 
ties are shown in Table 2 and are used in the combina- 
tions. 

Due to the additional event selection requirements 
used in the Id-analysis to optimise the expected un- 
certainty described in Section 5, for both channels 
the 2d-analysis has the smaller statistical uncertainty, 
despite the better top quark mass resolution of the 
Id-analysis. Both analyses are limited by the system- 
atic uncertainties, which have different relative contri- 
butions per source but are comparable in total size, 
i.e. the difference in total uncertainty between the most 
precise and the least precise of the four measurements 
is only 16%. 

The four individual results are all based on data 
from the first part of the 2011 data taking period. The 
e+jets and /i+jets channel analyses exploit exclusive 
event selections and consequently are statistically un- 
correlated within a given analysis. In contrast, for each 
lepton channel the data samples partly overlap, see Sec- 
tion 4. However, because the selection of the jet triplet 
and the construction of the estimator of mt op are differ- 
ent, the two analyses are less correlated than the about 
50% that would be expected from the overlap of events. 

The statistical correlation of the two results for each 
of the lepton channels is evaluated using the Monte 
Carlo method suggested in Ref . [39] , exploiting the large 
Monte Carlo signal samples. For all four measurements 
(two channels and two analyses) , five hundred indepen- 
dent pseudo-experiments are performed, ensuring that 
for every single pseudo-experiment the identical events 
are input to all measurements. The precision of the 
determined statistical correlations depends purely on 
the number of pseudo-experiments performed, and in 
particular, it is independent of the uncertainty of the 
measured mtop per pseudo-experiment. In this analysis, 
the precision amounts to approximately 4% absolute, 
i.e. this estimate is sufficiently precise that its impact on 
the uncertainty on mt op , given the low sensitivity of the 
combined results of mt op to the statistical correlation, 
is negligible. For the Id-analysis, the signal is comprised 
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of tt and single top quark production, whereas for the 
2d-analysis the single top quark production process is 
included in the background, see Table 1. Consequently, 
the MC@NLO samples generated at m top = 172.5 GeV 
for both processes are used appropriately for each anal- 
ysis in determining the statistical correlations. The sta- 
tistical correlation between the results of the two anal- 
yses is 0.15 (0.16) in the e+jets (/x+jets) channels, re- 
spectively. Given these correlations, the two measure- 
ments for each lepton channel are statistically consis- 
tent for both lepton flavours. 

The combinations of results are performed for the 
individual measurements and their uncertainties listed 
in Table 2 and using the formalism described in Refs. [39, 
40]. The statistical correlations described above are used. 
The correlations of systematic uncertainties assumed in 
the combinations fall into three classes. For the uncer- 
tainty in question the measurements are either consid- 
ered uncorrelated p = 0, fully correlated between anal- 



yses and lepton channels p — 1, or fully correlated be- 
tween analyses, but uncorrelated between lepton chan- 
nels denoted with p = (1). A correlation of p — is 
used for the sources method calibration and jet energy 
scale factor, which are of purely statistical nature. The 
sources with p = 1 are listed in Table 2. Finally, the 
sources with p = (1) are QCD background normalisa- 
tion and shape that are based on independent lepton 
fake rates in each lepton channel. 

Combining the results for the two lepton channels 
separately for each analysis gives the following results 
(note that these two analyses are correlated as described 
above) : 

m t0 p = 174.4 ± 0.9 sta t ± 2.5 syst GeV (Id-analysis), 
m t0 p = 174.5 ± 0.6 sta t ± 2.3 syst GeV (2d-analysis). 

For the Id-analysis the /i+jets channel is more precise, 
and consequently carries a larger weight in the com- 
bination, whereas for the 2d-analysis this is reversed. 



Measurement of the Top Quark Mass with the Template Method using ATLAS Data 



17 



u. 1.02 
CO 

—3 

1.01 

1 
0.99 
0.98 



0.97 



; 


— 1a 




ATLAS 


; 


\ 


— 2<5 




JSF = 0.985 ±0.008 

m, op = 1 74.3 ±1.0 MF GeV 


; 


- 


e + jets 


5 


s 


'» 


- 


rj 


L dt = 1.04 lb" 1 




, I , 


. , , 1 , , , . I . , . , I , . 


-_ 



171 172 173 174 175 176 177 178 179 

™,o P [GeV] 
(a) e+jets channel 



u. 1.02 
CO 



1.01 : 

i- 

0.99; 

0.98 : 
0.97 



1o 
2a 



JSF = 0.986 ±0.006 

m, op = 1 75.0 ±0.9 sla „ JSF GeV 



VRS\ 



H + lets 
L dt = 1.04 lb" 1 



171 172 173 174 175 176 177 178 179 

™,op[GeV] 
(b) /i+jets channel 
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However, for both analyses, the improvement on the 
more precise estimate by the combination is moderate, 
i.e. a few percent, see Table 2. 

The pairwise correlation of the four individual re- 
sults range from 0.63 to 0.77, with the smallest cor- 
relation between the results from the different lepton 
channels of the different analyses, and the largest cor- 
relation between the ones from the two lepton channels 
within an individual analysis. The combination of all 
four measurements of mtop yields statistical and sys- 
tematic uncertainties on the top quark mass of 0.6 GeV 
and 2.3 GeV, respectively. Presently this combination 
does not improve the precision of the measured top 
quark mass from the 2d-analysis, which has the better 
expected total uncertainty. Therefore, the result from 
the 2d-analysis is presented as the final result. The two 
analyses will differently profit from progress on the in- 
dividual systematic uncertainties, which can be fully 
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exploited by the method to estimate the statistical cor- 
relation of different estimators of mtop obtained in the 
same data sample together with the outlined combina- 
tion procedure. The results are summarised in Figure 10 
and compared to selected measurements from the Teva- 
tron experiments. 



9 Summary and conclusion 

The top quark mass has been measured directly via two 
implementations of the template method in the e+jets 
and /i+jets decay channels, based on proton-proton col- 
lision data from 2011 corresponding to an integrated lu- 
minosity of about 1.04 fb _1 . The two analyses mitigate 
the impact of the three largest systematic uncertain- 
ties on the measured TO t0 p with different methods. The 
e+jets and /i+jets channels, and both analyses, lead to 
consistent results within their correlated uncertainties. 
A combined Id-analysis and 2d-analysis result does 
not currently improve the precision of the measured 
top quark mass from the 2d-analysis and hence the 
2d-analysis result is presented as the final result: 



m to p 



174.5 ± 0.6 Bt at ± 2.3 syst GeV 



This result is statistically as precise as the mt op mea- 
surement obtained in the Tevatron combination, but 
the total uncertainty, dominated by systematic effects, 
is still significantly larger. In this result, the three most 
important sources of systematic uncertainty are from 
the relative 6-jet to light jet energy scale, the modelling 
of initial and final state QCD radiation, and the light 
quark jet energy scale. These sources account for about 
85% of the total systematic uncertainty. 
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